Recovery temperature for nonclassical energy transfer in atom-surface scattering 
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Nonperturbative expressions are derived for the angular 
resolved energy transfer spectra in the quantum regime of 
multiphonon scattering of inert gas atoms from surfaces. Ap- 
plication to He atom scattering from a prototype heatbath 
Xe/Cu(lll) shows good agreement with experiments. This 
enables a full quantum calculation of the total energy trans- 
fer \i and therefrom the much debated recovery or equilibrium 
temperature T r characteristic of zero energy transfer in gas- 
surface collisions in the free molecular flow regime. Classical 
universal character of \i and T r is refuted. 

PACS numbers: 68.35. Ja, 34.50.Dy, 63.22.+m, 47.45.Nd 



Despite the great progress recently made in under- 
standing the dynamics of gas-surface collisions |jj little 
effort has gone into applying this knowledge to problems 
of general relevance 0. Because of their fundamental 
importance for a wide range of the flow phenomena 
there is a need to be able to predict the magnitudes of 
the heat transfer fi [Q and the recovery temperature T r 
at which zero energy transfer to the surface occurs in 
the free molecular flow regime |5j , but the classically cal- 
culated (i and T r generally fail to reproduce the exper- 
imental data 0. However, the information on single- 
and multi-phonon excitations recently accumulated from 
He atom scattering (HAS) from surfaces |rj-|ic|l combined 
with novel theoretical developments now opens up the 
possibility of calculating the heat transfer on a micro- 
scopic level and in the fully quantum scattering regime. 
In this Letter a new approach based on multiphonon scat- 
tering theory [pd|-^3| is developed to predict the energy 
transfer in interactions of He atom beams with Xe mono- 
layers on Cu(lll) which represent an ideal prototype 
heatbath encompassing both the dispersive and nondis- 
persive surface phonons |i"o[ ] . The calculations reproduce 
the HAS data remarkably well and thus enable a reliable 
prediction of the recovery temperature which is found to 
differ markedly from the results of classical accommoda- 
tion theories currently in use 

The total energy transfer [i which enters the heat trans- 
fer and accommodation coefficients fEH is evaluated from 



£-^ki,T a (e)rfe- 



(1) 



perature T s . The final atom state can be either a con- 
tinuum | c) or a bound state | b) of the static atom- 
surface potential U(r) [p^| - [l6| . However, there are two 
major difficulties in applying Eq. (|l|) to atom-surface 
scattering. First, N^ it T B (£) and thereby /i(ki,T s ) are not 
directly accessible in typical HAS time-of-flight (TOF) 
measurements from which most of the data are available 
at present. These experiments yield the energy and an- 
gular resolved quantities usually only for fixed total scat- 
tering angle Osd = #j + Of, where Oi (Of) is the initial 
(final) scattering angle. The TOF spectrum is directly 
proportional to the c — > c component of the full energy 
and parallel momentum resolved scattering distribution 
Nk lt T s (e, AK) jb|] which in turn is related to the angular 
resolved energy transfer 



(j, r (ki,T s ,0f) 



fN™(e,AK(e))de ' 



(2) 



where N^ u t.(^) is the scattering spectrum which gives 
the probability density for an atom with initial momen- 
tum Kki to exchange energy e with a surface at the tem- 



As fi r (ki,T s ,0 f) can be computed from both the theo- 
retical N£Tt ( £ j AK) and the experimental TOF-spectra 
a direct comparison of the two results enables the verifi- 
cation of model calculations. The second difficulty arises 
in calculations of the reliable multiphonon N^ u T a (e), ob- 
tained from Nk it T,(s, AK) by integration over AK, in 
the regime in which it is important to treat the dynam- 
ics of both the projectile and surface vibrations quan- 
tum mechanically. However, the recent progress in in- 
terpreting the multiphonon HAS from monolayers of no- 
ble gas atoms (Ar, Kr, Xe) adsorbed on metals l^-|ic|l 
makes it now possible to exactly calculate N^ ;y T s (z, AK) 
and thereby assess the total energy transfer in the stud- 
ied collisions. These adlayers sustain low-energy longitu- 
dinal (L) and shear-horizontal (SH) in-plane polarized 
modes with acoustic-like dispersion and nondispersive 
vertically polarized Einstein- like modes (S). The presence 
of nondispersive modes gives rise to special interference 
effects in iVk, T,( e j AK) which requires its calculation in 
a closed, nonperturbative form. Such a solution is pre- 
sented below and applied to the study of energy trans- 
fer in HAS from (\/3 x ^/3)R30° monolayers of Xe on 
Cu(IIl) which exhibit very weak diffraction and hence 
can be considered as flat in energy exchange processes 
pot . This suppresses selective adsorption assisted energy 
transfer Jj"5[ and greatly simplifies calculations and com- 
parisons of N£~£ (e, AK) with the TOF spectra. 

In the scattering regime in which uncorrelated phonon 
processes are dominant the angular resolved scattering 
spectrum has the following unitary form ]Il|jL2| : 



1 



7V kiiTs (e,AK) 



drd 2 TL 



i[er— R(AK)R]/ft 



, (27T^) 3 

x exp[2W^(r,R) - 2W] 



(3) 



present quantum scattering regime the recoil is large (c.f. 
Figs. 10 and 11 in Ref. pi ) making ^(ki,T s ) strongly 
T s -dependent. 



Here 2W(t,H) is the scattering function Q, 2PF = 
2W(t — 0, R = 0) is the exponent of the Debye- Waller 
factor (DWF) |l7| and r and R are auxiliary integra- 
tion variables used to project the states with e and AK 
from the integral on the RHS of Eq. (|^). In the range 
of validity of Eq. (|^) the energy and parallel momen- 
tum are conserved in each phonon exchange process and 
e and AK are constrained to satisfy the total energy and 
parallel momentum conservation in the collision. Using 
expression (||) it is possible to write Eq. (Q) as: 



M(ki,T s ) = i^2W[r = 0,R = 0), 



(4) 



which can be readily calculated once 2W(t, R) is estab- 
lished. In the following the projectile-phonon coupling is 
assumed linear in the adsorbate displacements since this 
gives the dominant multiphonon contribution observed 
in HAS ftf]. This yields @: 

2W(r,R)= [l^£w G (+)| 2 [^Q,0 + l] 



x e 



^q 3 -t-(Q+g)R] 



V; 



Ki,Q+G 

k z ,k z i,j 



(-) I 2 n(uQj)e 



;[w qi T-(Q+G)R] 



(5) 



where Q, j and ujQj denote the parallel wave-vector, 
branch index and frequency of a normal phonon mode, 
respectively, and G is the adlayer reciprocal lattice vec- 
tor, k = (K, k z ) where k z is the quantum number de- 
scribing projectile perpendicular motion in c- or 6-states 
of U(r). n(uiQj) is the Bose-Einstein distribution, and 
V^'j|7 4 "j (±) denote one-phonon emission (+) and ab- 
sorption (— ) on-shell scattering matrix elements normal- 
ized to particle current in the z-direction ]lO| , |T2] | . Substi- 
tution of expression (||) into (^) gives 

//(ki,T s ) = Mo(ki) +/i rec (ki,T s ), (6) 

in which the temperature independent part is 

/ i (k i )= *»<v i i 2 . (7) 



FIG. 1. Comparison of the temperature dependence of 
the angular resolved energy transfer ^i r (ki, T s , Of) calculated 
from He— >Xe/Cu(lll) TOF spectra for four experimental Ei 
and fixed scattering geometry (open symbols), and from the 
present model (full lines) . The inset shows a comparison of the 
measured and calculated multiphonon scattering spectrum for 
Ei = 45.11 meV, 9 l = 50°, T s = 58.2 K. 

The nonperturbative 
multiphonon solution for JV^^^AK), which enables 
explicit calculation of expression (||) and thus a direct 
comparison with the HAS data, is obtained by separat- 
ing the most strongly coupling Einstein branch j = S of 
frequency lus out of the sum in Eq. (|^), giving: 

exp[2IF(r, R) - 2W] = (r, R)A^t s (t, R), (9) 

where dis denotes the remaining dispersive modes. Em- 
ploying trigonometric identities to N^j- (t, R) yields: 



^ — HUJST 



= — oo 



(10) 



where 2Ws — 2Ws(t = 0,R = 0) is the corresponding 
Debye- Waller exponent and 



and the T s -dependence is determined by the recoil term P[ (R) = 



/i rec (ki,T s ) = ^2 huQj 

Q,G,j,k, 



V, 



Ki,Q+G 

k z ,k z i,j 



( + ) I 5 




(w s ) + l]V|(R,+) 



n(w s )V|(R,- 



- 1 <:£!f (-) i 2 ] ^ 

which vanishes in the recoilless trajectory approxima- 
tion (TA) for the projectile motion |19[]. However, the 
TA may fail even for heavier atoms p0|,|l6| and in the 



v /4n( Ws )[n(o; s ) + l]V|(R, 



)Vi(R- 



where Vl(R,±) = £ Q , G , fe , | V££g*(±) 
eT i["sr-(Q+G)R] and /; is tlie mo dified Bessel function 
of the first kind. Hence, the separated Einstein phonon 
component of the spectrum in Eq. (0) takes the form 



(11) 



2 



J2 Ni{AK)S{e-lfius), (12) 



l=—oo 



where A^AK) = / d 2 Ke- i( - AK ^ R Pi(Ti)/(2w) 2 is the in- 
tensity of the Z-th Einstein multiphonon loss (/ > 0) 
or gain (I < 0) peak. The elastic intensity N®™ (e = 
0, AK 7^ 0) is non-vanishing for T s > because multi- 
ple exchange of nondispersive phonons may give rise to 
finite momentum transfer without a net energy trans- 
fer. However, the spectral intensity of the specular elas- 
tic peak is e~ 2Ws S(AK) because P (R oo) -> 1. On 
the other hand, the elastic peak in iVjk, t s (e) , obtained 
from Eq. (12) by integrating over all AK, is weighed 
by e~ 2Ws P a (R = 0) but with P a (R = 0) ^ 1 due to the 
same multiquantum exchange effect. 

A similar procedure yields for N^ S T (r, R) in Eq. (|j): 



two highest energy curves in Fig. || and confirmed by ad- 
ditional calculations at high Ei, manifests itself only for 
fixed 9i because the ^-dependence of three-dimensional 
scattering matrix elements is not contained solely in the 
factorizable scaling factor E Z i Also, extension of 

the classical Baule expression pertinent to energy transfer 
in the cubes model (AK = 0) to quantum surface scat- 
tering is generally inadequate, as is demonstrated in 
comparison with our results for T s = 0. Since U(z) with 
the well depth of 6.6 meV supports three bound states 
we show in the inset the low-E^ dependence of the recov- 
ery temperature (for which /i(kj,T r ) = 0) calculated for 
the present phonon heatbath for fixed ftj, for He coupling 
either only to S- or to S-, L- and SH-phonon modes. The 
small difference indicates that the major contribution is 
from strong He atom coupling to vertically polarized S- 



N^ S T (t,R) = e~ 2W/ ^exp 



Q,G,fe' 8 ,&S 



x cos \ i In , 



[n( UQj ) + 1] | V, 



\ «("<w) I V££+ G R I 



Expressions (|Io|)-(|l3|) are exact within the validity of (||) 
fl2"f and include the combined effect of recoil and temper- 
ature on the multiphonon scattering spectra. 

The calculations of expressions (||)-([l3|) were carried 
out by modeling U(r) = U(z), the projectile-phonon 
coupling and the dynamical matrix of the Xc monolayer 
and 40 layer thick Cu(lll) slab as in Refs. §|uj. The 
calculated N£~tp (s, AK) and /x r (ki, T s , Of) in the mul- 
tiphonon scattering regime were tested by directly com- 
paring with experiment. Figure |l| shows a comparison of 
the angular resolved energy transfer obtained from Eq. 
(||) without invoking any fit parameters, with the values 
calculated by integrating the TOF spectra for fixed 6sd 
P.p"o| . The inset shows a comparison of the experimental 
and calculated scattering spectrum (Q) for one particular 
set of the scattering parameters. Thus confirmed consis- 
tency of the theoretical with experimental results enables 
consistent calculation of /i(ks,T s ). 

Figure [2] shows the temperature dependence of total 
heat transfer in HAS from Xe/Cu(lll) normalized to ver- 
tical component of the projectile incident energy, E Z i = 
Ei cos 2 9{. The T s -dependence of /i rec (kj, T s ) hinders en- 
ergy transfer to phonons and causes negative slopes in the 
plots. This arises from larger phase space for projectile 
c — > c transitions into final states with E z > E Z i. There 
it may give rise to negative /Lt(ki,T s ) (e.g. for Ei = 2.4 
meV and 9i — 50° at T s > 62 K) and hence to heating of 
the scattered beam. In the classical theory this effect is 
independent of the accommodation coefficient and hence 
of 9i ||. Here, the universal behaviour of |ii(kj,T s ) for 
higher Ei , as exemplified by the near coincidence of the 



[w Qj T - (Q + G)R] 



Ki,Q+G 



(+) I 2 



(13) 



modes E0 |. Rapid variations in T r are caused by kine- 
matic focusing in S-phonon assisted c — > b transitions for 
large parallel momentum transfer. 



FIG. 2. Temperature dependence of the total energy 
transfer [i in He^Xe/Cu(lll) collisions normalized to E z i 
for 9i = 50° and Ei = 80 meV (full curve), 60 meV 
(long-dashed curve), 20 meV (short-dashed curve) and 2.4 
meV (dash-dotted curve). The corresponding normalized val- 
ues of [i at T 3 = obtained from the Baule formula corrected 
for the well depth are denoted by inverted triangle, circle, dia- 
mond and triangle, respectively. The scattered beam is heated 
on the average if jj, < 0. Inset: relative contributions of the 
phonon modes to the recovery temperature T r (see text). 
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For monoatomic gases the standard classical accom- 
modation theory gives T r — Ei/2ks and to a good ap- 
proximation Ei ~ 5fcsTo/2 where To is the stagnation 
temperature of beam gas prior to expansion in the noz- 
zle This yields the recovery factor T t /Tq ~ 1.25 but 
deviations from this universal behaviour have been ob- 
served in wind tunnel molecular beam experiments Q 
and their explanation was proposed in terms of heuristi- 
cally modified classical expressions || . The present quan- 
tum theory enables essential progress beyond the classical 
results by allowing the parallel momentum exchange with 
phonons, multiphonon interference and quantum recoil of 
the projectile. Their interplay gives the recovery factor 
as a function of Ei (or To) and 9i which for our proto- 
type heatbath is shown in Fig. ||. Quite generally, T r is 
largest for normal incidence and only at higher Ei (i.e. 
To) quantum results may approach the classical limit so 
far observed only for rough technical surfaces Q|. Large 
deviations from the classical limit at low Ei are due to 
the quantum regime which allows larger AK and transi- 
tions affected by the bound states of He-surface poten- 
tial. Qualitatively these findings are not system specific 
as the present theory is quite general and can be readily 
extended to calculations of the heat transfer in collisions 
of He [^l],^2| or heavier rare gas atoms |ll|^3| with clean 
surfaces in a wide range of scattering conditions. 



The work in Zagreb has been supported in part by the 
NSF grant JF-133. 



FIG. 3. Recovery factor T r /To for a prototype heatbath 
sustaining phonon modes typical of Xe/Cu(lll) system plot- 
ted as a function of the incident energy Ei or stagnation tem- 
perature To of the gas for three representative incident angles 
9i. Dashed-dotted line is the classical result T r /To — 1.25. 

In summary, we have shown that in inelastic gas- 
surface scattering under the conditions of free molecular 
flow the combination of quantum and temperature effects 
gives rise to a violation of the universality of the energy 
transfer and recovery factor predicted by the classical ac- 
commodation theory. 
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